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Abstract. Multisite phosphorylation networks are encountered in many intracellular processes 
like signal transduction, cell-cycle control or nuclear signal integration. In this contribution net- 
works describing the phosphorylation and dephosphorylation of a protein at n sites in a sequential 
distributive mechanism are considered. Multistationarity (i.e. the existence of at least two pos- 
itive steady state solutions of the associated polynomial dynamical system) has been analyzed 
and established in several contributions. It is, for example, known that there exist values for the 

C^ ' rate constants where multistationarity occurs. However, nothing else is known about these rate 

^~H , constants. 

^3 ' Here we present a sign condition that is necessary and sufficient for multistationarity in n- 

site sequential, distributive phosphorylation. We express this sign condition in terms of linear 
^^. systems and show that solutions of these systems define rate constants where multistationarity is 

C^ ' possible. We then present, for n > 2, a collection oi feasible linear systems and hence give a new 

and independent proof that multistationarity is possible for n > 2. Moreover, our results allow 
to explicitly obtain values for the rate constants where multistationarity is possible. Hence we 

C^^ ' believe that, for the first time, a systematic exploration of the region in parameter space where 

multistationarity occurs has become possible. 

Keywords: sequential distributed phosphorylation; mass-action kinetics; multistationarity; sign 
condition; rate constants 

d 

1. Introduction 

CN \ Protein phosphorylation and dephosphorylation is frequently encountered in procaryotic and in 

^ ' eukaryotic cells. Thereby the majority of proteins is phosphorylated at more than one phosphoryla- 

vQ . tion site (with most procaryotic proteins having n <7 phosphorylation sites and eukaryotic proteins 

^3 ! generally having a large number of phosphorylation sites - certain eukaryotic proteins have up to 

^O ■ n = 150 sites |26j). Multisite phosphorylation networks are consequently encountered in many in- 

■rij" , tracellular processes, including, among others, the following: (i) signal transduction networks from 

^^ ' the family of Mitogen Activated Kinases (MAPK cascades), where the core module consists of one 

protein phosphorylated at a single site and two proteins each phosphorylated at two sites (see e.g. 

|16 tl23l[2Hl31| ): (ii) check-points in cell-cycle control, for example the control of the Gl/S transition 

in budding yeast involving the protein Sid with n = 9 (see e.g. [T1[21[2H]) and (iii) nuclear signal 

^\f . integration, for example, in mammalian cells involving members of the A/^F^IT transcription factor 

^; family with n= 13 (see e.g. mill OdHl [27]). 

" ■ ■ In this contribution we focus on networks describing the phosphorylation and dephosphorylation of 

a single protein A at n sites in a sequential distributive mechanism (cf. Figure [1] and, for example, 

|13|). Such a mechanism is able to cover the essential aspects of proteins like the aforementioned 

Sid or the members of the NFAT family |2H I22| . 

Due to the biological significance of sequential distributive phosphorylation, mathematical models 
have been studied in a variety of publications. Several publications deal with double phosphoryla- 
tion (i.e. n = 2) and establish multistationarity [ll[5l[6l[TT] (i.e. the existence of at least two positive 
steady state solutions, cf. Definition I3.5p and bistability |19| (i.e. the existence of two (asymptoti- 
cally) stable steady state solutions). Networks with n phosphorylation sites have been studied, for 
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Figure 1. Network describing the sequential distributive phosphorylation and de- 
phosphorylation of protein A at n-sites by a kinase Ei a phosphatase E2. The 
phosphorylated forms of A are denoted by the subscript nP (denoting the number 
of phosphorylated sites) . 



example, in [21|, 122]. where bistability is demonstrated numerically and in [12j, where it is argued 
that even though bistability is in principle possible, it is unlikely to occur in-vivo. 
The set of positive steady states of networks with n phosphorylation sites has been studied alge- 
braically in [25] and in [20] . In [25] the set of positive steady states of the more general class of so 
called post translational modification networks (with n-site sequential distributive phosphorylation 
as a special case) is analyzed. Such networks admit a rational parameterization of the set of pos- 
itive steady states. And in [20] it is shown that for n-site sequential distributive phosphorylation 
this set is defined by a binomial ideal. Hence n-site sequential distributive phosphorylation is an 
instance of a Chemical Reaction System with toric steady states. In [20] a condition for multista- 
tionarity in these reaction systems is given in terms of a sign condition on elements of two linear 
subspaces. Based on the results presented in ^0], this condition can be restated in terms of linear 
systems: to decide the sign condition, a large set of linear systems has to be examined and multi- 
stationarity is possible, if and only if at least one of these is feasible (cf. [2^ Theorem 5.5] and |31 
Lemma 2 & Theorem 2]). 

The number of steady states in general n-site sequential distributive phosphorylation has been 
studied in |30| . There upper and lower bounds have been established: there exist parameter values 
such that, for n even (odd), there are at least n+1 (n) steady states. And for all (positive) parameter 
values there are at most 2 n — 1 steady states. Hence multistationarity has been established for 
sequential distributive phosphorylation in |30j. However, no information other than existence is 
given about the parameter values where multistationarity is possible. In this contribution we address 
this question: we present for every n > 2 a collection of feasible linear systems and show that 
every solution of one of those systems defines parameter values where multistationarity is possible 
(together with two positive steady states as witness). Thus these results not only constitute a new 
and independent proof that multistationarity is possible for n > 2, they additionally enable, for 
the first time, a systematic exploration of the region in parameter space where multistationarity 
is possible. Due to the ubiquity of multisite phosphorylation this may be of potential interest to 
researchers working in many fields of (quantitative) biology. For the purpose of finding parameter 
values in biologically meaningful ranges we will start with this exploration in |15| . 
We arrive at our main result on the basis of our previous work from [S", ^ and |20| : as in the 
aforementioned references we derive a sign condition that is necessary and sufficient for the existence 
of multistationarity (Theorem 14.21 fc l4.6p . In contrast to |3l|2^ our condition is formulated in terms 
of the sign patterns generated by two linear subspaces (given as the image of two matrices defined 
in eqns. ([7e|) & ()15cp ). We then exploit the fact that every sign pattern uniquely defines a linear 
system to restate the sign condition in terms of these linear systems: multistationarity is possible, 
if and only if at least one out of ^ (3(3"+3) _ 1) or 2*-^""'"^^ linear systems is feasible depending on 
whether some or all components of the two witness steady state vectors have to differ to qualify 



MULTISTATIONARITY IN SEQUENTIAL DISTRIBUTED MULTISITE PHOSPHORYLATION NETWORKS 3 

for multistationarity, cf. Definition 13.51 and Remark 14.71 For the latter case, we present formulae to 
construct all sign patterns defining feasible linear systems (cf. Theorem 14. 9p and show that every 
solution of one of these linear systems defines parameter values where multistationarity occurs with 
two positive steady states differing in every component. We show that there are 2 (n — l)(n + 2) 
such sign patterns (cf . Proposition I4.12p . 

This paper is organized as follows: Section [5] introduces some of the basic notations. In Section [3] 
a mass action network is derived from Fig. [1] and the associated dynamical system is defined along 
with a formal definition of multistationarity. In Section |3]the main results address multistationarity 
in terms of sign patterns defining feasible linear systems (cf . Theorem 14.21 & I4.9p . The proofs of two 
of these results have been relegated to separate Sections [5] & [6] for easier reading. 



2. Notation 

We use the symbol M^ to denote Euclidian m-space, the symbol iR>o *° denote the nonnegative 
orthant and M'^q to denote the interior of the nonnegative orthant. Vectors are considered as column 
vectors and, for convenience, usually displayed as row vectors using ^ to denote the transpose. For 
example, x € JR"^ will usually be displayed as (xi, . . . , Xm) ■ 
For positive vectors x G -ff?>o we use the shorthand notation Inx to denote 

Inx := (Inxi, . . . , InXm) ■ (la) 

Similarly, for x € M"^, we use e^ to denote 

e^:=(e^i, ...,e^™)^ (lb) 

and, for x G JR^ with Xi ^ 0, i = 1, . . ., m 

x-^ :=(-,..., -)^ (Ic) 

Finally, x^ with x, y € JR>q will be defined by 

m 

xy:=J{xf. (Id) 

3. Mathematical models of n-siTE sequential distributive phosphorylation 

Here we first introduce the notation used to describe dynamical systems defined by mass action 
networks by means of the example network depicted in Fig. [2l Then we discuss the mass action 
network derived from Fig{T]and present the dynamical system defined by this mass action network. 

3.1. Dynamical systems defined by a single phosphorylation network. The network de- 
picted in Fig. [2] describes the phosphorylation of a protein A at a single site by the kinase Ei and 
the dephosphorylation of mono-phosphorylated protein Ap by the phosphatase E2 (cf. Fig. [2aj) . Fig- 
ure [2b] contains a mass action network describing this process, where chemical reactions are denoted 
by arrows pointing from the species that are consumed in a reaction (the educts) to the species that 
are produced (the products) and fcj (/j) denote rate constants. The network consists of 6 species, the 
kinase Ei (with concentration variable xi), the protein A (^2), the phosphatase E2 (a^s), the enzyme 
substrate complexes AEi (X4) and ApE2 (xg) and the mono-phosphorylated Ap (xs). The network 
is composed of 6 reactions, two describing the reversible binding of protein and kinase Ei+A^=^AEi, 



KATHARINA HOLSTEIN, DIETRICH FLOCKERZI AND CARSTEN CONRADI 
P 




(a) 1-site phosphorylation 



xi X2 ki X4 fcs xi X5 

Ei + A . AEi >Ei + Ap 



ko 



h 



X3 X5 ll X6 t3 X3 X2 

E2 + Ap ^=^ ApE2 ^E2 + A 



(b) Mass action network 



Figure 2. Phosphorylation of protein A at a single phosphorylation site: process 
diagram (a) and mass action network (b). 



one the irreversible release of mono-phosphorylated protein and kinase AEi^Ei + Ap, two the re- 
versible binding of mono-phosphorylated protein and phosphatase E2 + Ap^ApE2 and one the 
irreversible release unphosphorylated protein and phosphatase ApE2^^E2 + A. 
To every reaction we associate a reaction vector describing the stoichiometry of that reaction. For 
example in the reaction 

El + A >AEi 

one unit of A {X2) and Ei [xi) are consumed and one unit of AEi (X4) is produced, which yields 
the vector 

(-1, -1, 0, 1, 0, of . 

All in all one obtains six reaction vectors, that we collect as column vectors of the stoichiometric 
matrix 

" -1 1 1 
-110 1 
0-111 
1-1-1 
1-110 
1-1-1 



5(1) 



where the superscript '^^ is used to denote the stoichiometric matrix of a 1-site phosphorylation 
network. Columns 1 and 2 of S^^> correspond to the reversible binding of A and Ei, column 3 to 
the irreversible release of Ap and Ei , columns 4 and 5 to the reversible binding of Ap and E2 and 
column 6 to the irreversible release of A and £'2- 

Using mass action kinetics the reaction rate Vi (i.e. the speed of reaction i) is proportional to the 
product of the concentrations of the educts, for example, the rate ri = fci xi X2 for the reaction 
Ei+A^ AEi . Similarly one obtains r2 = k2 X4 for AEi -^ Ei + A, r^ = k^ X4 for AEi -^ Ei + Ap, 
r4 = h X3 X5 for E2 + Ap —^ ApE2, r^ = h xq for ApE2 —^E2 + Ap and r^ = I3 xq for ApE2 -^ E2 + A. 



iT 



(fci, ..., fcs, /i, ..., hy) 



We collect ri, . . . , rg in a vector (with x = (xi, . . . , xq 

r(i) {k, x) = [ki xi X2, k2 X4, k^ X4, h X3 X5, h xq, k X(,f 
and obtain the dynamical system (where we suppress the time dependence of Xj): 



(2) 
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with 

Xl = {k2 + k^) X4 - kiXiX2, 
±2 = -kiXiX2 + k2X4 + /3X6, 
^3 = (^2 + ^3)a;6 -^1X3X5, 

±4 = - (/c2 + fca) X4 + kiXiX2, 
X5 = /C3X4 - /1X3X5 + /2X6, 

X6 = - ih + ^3) X6 + /1X3X5. 

Note that the vector A^' {k, x) is fully characterized by the vector k and the exponents of the 
monomials. We collect those in the rate exponent matrix y'^^': 
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3;(i) 



Fig.[2]does not consider protein synthesis and degradation, the total amount of protein A is therefore 
constant, that is, one has the conservation relation 

X2 + X4 + X5 +X6 = Ci, 

where ci is a constant denoting the total concentration of A. The name conservation relation stems 
form the fact that the above sum of concentrations is constant along solutions of ([2]): from ^ it is 
easy to see that ±2 (t) + X4 (t) + X5 (t) + xq (t) = and hence X2 (t) + X4 (t) + X5 (t) + xq (t) = const. 
Likewise one has 

xi + X4 = C2 and X3 + xg = C3, 

with constants C2 and C3 denoting the total concentration of Ei and E2, respectively. Conservation 
relations are defined by elements of the left kernel of the matrix S^^' . One obtains, for example, the 
following full row-rank matrix 



1 








1 














1 








1 





1 





1 


1 


1 



^(1) 



with Z''^' S^^' = where the rows form a basis of the left kernel of S^^' . 

3.2. The mass action network derived from Figure [Tl The mass action network derived 
from Fig. [1] (with n an arbitrary but fixed positive number) has been described in |20| I30j . Here 
we use a similar mathematical description: the network consists of the following chemical species: 
the protein (substrate) A together with n phosphoforms Ap, . . . , Anp; the kinase Ei together 
with n kinase-substrate complexes AEi, . . . , An-ip Ei and the phosphatase E2 together with n 
phosphatase-substrate complexes Ap E2, . . . , AnP £'2- Hence there is a total of 3n + 3 species. To 
each species, a variable Xi denoting its concentration is assigned as depicted in Table [TJ We collect 
all variables in a (3n + 3)-dimensional vector 



(xi, 



X3n+3j 



As it will turn out, the chosen labeling entails a simple block structure for the matrices associated 
to the dynamical system (jlOap of the network in Fig. [H cf., for example, the block structure ()13bp 
for the generators of the nonnegative cone in the kernel of the stoichiometric matrix. 
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Phos. # 


Species 


Var. 





El 
A 
E2 


Xi 
X2 
X3 


1 


A El 

Ap 

ApE2 


X4 
X5 

xe 








i 


Ai-ip El 
AiP 

ApE2 


Xl+3i 
X2+3i 
X3+3i 








n 


An-lP El 

AnP 

AnpE2 


Xl+3n 
X2+3n 
X3+3n 



Table 1 . Assignment of variables to species 



Assuming a distributive mechanism, a single phosphorylation occurs with each encounter of sub- 
strate and kinase [13] and n phosphorylations therefore require n encounters of substrate and kinase 
of the form 

El + Ai^ip ^^=^ Ai^ip El >Ei + AiP, i = 1, . . . , n , 

where Aqp = A and Aip = Ap. Similarly, n dephosphorylations following a distributive mechanism 
require n encounters of substrate and phosphatase of the form: 



E2 + AiP ^=^ AiP E2 >E2 + Ai_ip, 



1,... ,n . 



Each phosphorylation and each dephosphorylation therefore consists of 3 reactions and consequently 
the network consists of 6n reactions. To each reaction we associate a rate constant. We use ki for 
phosphorylation and li for dephosphorylation reactions and obtain the following reaction network: 



k3i-2 hi 
El + Ai^ip V Ai^ip El >Ei + AiP, i = l, ...,n 



E2 + AiP 



hi-i 

hi-2 



(3) 



3i 



=i^ AiP E2 >E2 + vlj-ip, 1 = 1, ... ,n 



3i-l 



Using this notation, /c3i_2 (^31-2) denotes the association constant, /c3j_i (^si-i) the dissociation 
constant and k3i (Zsj) the catalytic constant of the i-th phosphorylation (dephosphorylation) step. 
We collect all rate constants in a vector k defined in the following way: 

Definition 3.1 (The vector of rate constants k). 

We collect the six rate constants associated to the i-th phosphorylation /dephosphorylation in net- 
work ^ in a sub-vector 



'^{i) '■— (^3i-2, fcsi-l, ksi, l3i-2, hi-l, hi) 

We then combine the sub-vectors to a 6n- dimensional column vector 

K :=co1(k(i), ..., K(„)) . 



(4a) 



(4b) 
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We conclude this subsection with a brief comment on the notation used to denote (sub-)vectors and 
matrices throughout this contribution. 

Remark 3.2 (Vector and matrix notation). 

In the foUowing sections the dimension of vectors is determined by the number of phosphorylation 
sites n according to the formula 3n + 3. (The vector x, for example, lives in Euclidian IB?^~^^, cf. 
Table [T|) . We will use the symbol Cj to denote elements of the standard basis of Euclidian vector 
spaces and use the superscript '*•* to distinguish basis vectors of vector spaces of different dimension 
3i + 3: 



,(^) 



, denotes elements of the standard basis of ]R 



3i+3 



(5) 



Likewise for matrices the superscript ^^' is used to indicate that the number of rows and/or columns 
depends on an integer i (see, for example, equations ([9]), ([7e|l . ()8dp or (|13bp below). 
Later on we will split vectors of length 3n + 3 into consecutive sub- vectors of length three and we 
will use the subscript u\ to denote the i-th sub-vector 



X(i) ■= (Xsi+l, X3i+2, Xsi+^y 



and we will use x = col {x(q\, . . . , Xfn)) to denote that x consists of n+ 1 such sub- vectors. Likewise 
we will split vectors of length 6n into sub- vectors 



2/(i) := {y6i~5, • • • , yeiY 



of length six and use y = col (^(i), • • • , ?/(«)) to denote that y consists of n such sub- vectors. □ 



3.3. The dynamical system defined by the mass action network ([3]). Starting with the 
stoichiometric matrix S^^' and the rate exponent matrix y^^> defined by the mass action network 
depicted in Fig. [2b] one can recursively construct matrices S^"'' and y^"^' for the network ^. Using 
the ordering of species and reactions introduced above in Table [1] one obtains S^'^' G jR(3™+3)x6n^ 
y{n) g j^(3n+3)x6n ^nd z(") by the following steps: 

(I) Concerning S^^': 



5(1) 



nil 

^21 



5(2) 



" 5(1) 


ni2 

"22 


. 03x6 


"21 



(6a) 



SU) 



5(i-i) 











3x60-1) 



"12 

3(j-2)x6 
?^22 
"21 



(6b) 
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for j = 3, ..., n with the following sub-matrices of dimension 3x6: 



nil 



^12 



^21 



^22 



For fixed n, the recursive formula ()6bp evaluates to 



1 


1 


1 











1 


1 











1 











-1 


1 


1 


1 


1 


1 






































-1 


1 


1 


1 


-1 


-1 

















1 


-1 


1 














1 


-1 


-1 




















1 


1 











1 





















(7a) 
(7b) 
(7c) 
(7d) 



g(n) .__ 



nil 

^21 


"-12 
"■22 
"21 


"12 
03x6 
"22 
"21 


"12 

02-3x6 

"22 
"21 




03x6 


03x6-2 


03x6-3 









3x6-(n- 


1) 





(II) Concerning 3^("-): 



"12 



0(n-2)-3x6 



"21 



p(l) ,p(l) p(l) p(l) p(l)+p(l) p(l) p(l) 
^1 ^ ^2 ^A ^A ^3 ^ ^5 ^6 ^6 



yU~i) 



3x6(j-l) 



(i) , Ji) Jj) Jj) Ji) , Ji) Ji) Jj) 



eY + e 



^31+1 ^37+1 ^3 + ^37+2 ^37 



3j-l ^3j+l ^3j + 



3j+2 "=^3j+3 '='3j+3 



for j = 2, ..., n. Using 

Yoii:):-- 
the recursive formula (l8bD evaluates to 



^1 T^ ^3i-l '^3J+1 ^3i+l ^3 ^ ^3i+2 ^3i+3 ^3J+3 



(7e) 



(8a) 
(8b) 

(8c) 



yin)^ .,-_ 



Yo{lf 06-1x3 


06-2x3 


06-3x3 








Yoi^r 




yo{3f 


06-nx3 












Yo{n-iy 










lo(ny^ 









(8d) 



(III) Concerning Z^""' representing a basis for the left kernel of S'*^""-': 
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^H 



1 


1 




1 


1 


1 




1 


1 


1 1 1 




1 1 1 












n— times 



(9) 



Definition 3.3 (The dynamical system defined by network ([3])). 

The vector of rate constants k from ^J^ o^^^d the columns yi of 3^'""' define a monomial function 



r(") (k, x) = diag (k) (x^S . . . , x^«")^ =: diag (k) $(") (x) 



(10a) 



(cf. Iild\) ) and, together with the stoichiometric matrix S'"-* from ^7^, the dynamical system 

X = 5(") r(") (k, x) (10b) 

defined by the network (0) of Fig. QJ 

Remark 3.4 (The monomial function r'-"')(K, x)). 

Observe that the factor <!>("'•' (x) in the monomial function (jlOap satisfies 



$(") 



for vectors /i € M^ and 



„y("'^M 



$(") (x-1) 



(11a) 



lib) 



$(n) (^x) 
for vectors x E iR™ with Xj 7^ 0, i = 1, . . . , m. 

From equation (|10bp follows that level sets {x € iR'^"'^^ |Z '-"'•* x = const.} are invariant under the 
flow of (|10bp as Z*^"^ x(i) = Z^""^ x(0) along solutions x(t) of (|10bp . This observation motivates the 
classical definition of multistationarity originating in chemical engineering (cf. Figure [3] and, for 
example, [31 Definition l&Remark 2] or [9l I10|): 

Definition 3.5 (Multistationarity). 

The system x = S^"^' r^""' {k, x) from (|10bp is said to exhibit multistationarity if and only if there 

exist a positive vector k G IR^q and at least two distinct positive vectors a, b (z IRyo^ with 

5(") r(") {k, a) = (12a) 

5(") r(") {k, b) = (12b) 

ZH a = z(") b. (12c) 

We conclude this section by a discussion of ker (S'"-*) and of the pointed polyhedral cone ker (S'"-*) fl 

7p3n+3. 

Lemma 3.6 (A nonnegative basis of ker (S"*-"-*)). 

Let n > 2 be given and recall the stoichiometric matrix S^^' as defined in ^7^. Let 

1 1 

1 

1 

1 1 

1 

1 



E 



(13a) 
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Figure 3. Multistationarity and bistability: While the equation S^"'^ r(")(K,x) = 
often has infinitely many steady state solutions, a solution x{t) of the differential 
equation x = S^"^' r^"''{K,x) 'sees' only those contained in the level set L^tQ\ := {x G 
j^3n+3|2'{") 3; = Z''"'' x{0)}. Multistationarity requires at least two positive steady 
state solutions in L^^m), bistability two 'stable' positive steady state solutions in 



L 



x(0)- 



and define the 6n x 3n matrix with E on the block diagonal (and zero otherwise): 



^W ., 



E 



E 



(13b) 



i3n+3 



Then E^"^' is a basis o/ker (S"*^"-*) and the columns of E^""' are generators o/ker (S"*-"^) n IR^q . 

Proof. First recall the building blocks of S^"'\ the 3x6 sub-matrices nn, 7221, ni2 and 7222- We 
show that im (£'*^")) = ker (S**^")) by proving the inclusions im [E^"'^) C ker (S*^")) and ker (5*^"^) C 
im (£;(")). 

(i) Concerning im (£'(")) C ker (S^")): 

Straightforward computations show that nuE = O3X3, ni2E = O3X3, n2iE = 63x3 and n22 E = 

03x3- ThusSWi?W=0(3n+3)x3n- 

(ii) Concerning ker (S"*^")) C im [E^^^): 

First verify that i? is a basis of ker (n2i) and ker (n22)- Pick any vector r] € ker (S"*-"-*) and split 
it in consecutive sub- vectors of length 6 as described in Remark 13.21 r] = col (??(i), • • • , Vin))- The 
vector r] satisfies S^""' 7/ = 0, that is (cf. eq. ([7e|) defining S'"') 

n 

"■11^(1) +^'^l2Vii) =0, 
j=2 

n21??(l) + ^22^/(2) = 0, 



n2l1]{n~l) + n22??(n) = 0, 
"'21??(n) = 0. 
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From the last equation follows that t//^) G ker(n2i) = im.{E). The preceding equations then imply 
r/(j) G ker(n2i) = im (E), for i = n — 1, . . . , 1. Hence r] £ im (£^'"^) and thus ker (S'")) C im (£^'"^) 
(recall that r/ E ker (S"*-"-*) was chosen arbitrarily). 

Finally, concerning the pointed polyhedral cone, observe that a vector 7] is called a generator of 
the cone ker (S^"^) n Si^Q^ , if and only if t] satisfies the following three conditions (cf. for example, 

my- 

r?>0 

5(") r? = 

And if 7/1 and r/2 are generators of ker (S"*-"^) PI IR^q^ , then the following implication holds: 

supp (771) C supp (7/2) =^ ??i = or r/i = ar/2, with a > 0, 

where for vectors r/ S JR'^"+^ the symbol supp (77) denotes its support (i.e. the set of indices of the 
nonzero elements). Observe that the column vectors of E satisfy these conditions and recall the 
block diagonal structure of E^"-' . Hence the columns of E^""' satisfy these conditions as well. D 



4. Results 

In this section we first present a condition for multistationarity (cf. Definition 13. 5p in terms of sign 
patterns, followed by a discussion of sign patterns satisfying this condition. 



4.1. A sign condition for multistationarity in the dynamical system ()10b|) . We now turn 
to multistationarity for dynamical systems defined by network ^, that is to the system defined 
in (llObp . It follows from Definition 13.51 that multistationarity requires, for a given vector k, the 
existence of two positive solutions a and b of the polynomial equations (|12ap and (|12bp that satisfy 
the linear condition (|12cp . In Theorem 14.21 below, we first turn to the polynomial equations. The 
linear constraint is taken into account afterwards in CoroUarv 14.51 In the following Theorem 14.21 we 
use matrices II'"' and M'-"-' to argue that, for system (|10bp . the existence of positive solutions of 
the polynomial equations is equivalent to the existence of solutions of a linear system. 

Definition 4.1 (Matrices H^") and M^")). 

Let 1 denote a vector of dimension 6 filled with the number 1. Define the 6x2 matrix 

Ho (i) : = [ (2 - i) 1 (i - 1) 1 ] , i = 1, . . . , n (14a) 



and the 6n x 2 matrix 

Further define the 3x3 matrices 



Hnfl) 



no(n) 



(14b) 



M(0,7z) := 



and 



M{i,n) :-- 



-1 —n +1 n 

1 n —n 

-1 -n + 2 n-1 



-i + 2 7-1 

1 n — i —n + i 
-7 + 2 7-1 



(15a) 



(15b) 



12 KATHARINA HOLSTEIN, DIETRICH FLOCKERZI AND CARSTEN CONRADI 

for i = 1 , . . . , n and the (3 n + 3) x 3 matrix 

' M{0,n) 
M(-.):= ''<!■"' 

M {n,n) 
In Section [5] below we prove the following result. 



(15c) 



Theorem 4.2 (Solutions of the polynomial equations). 

Consider the n-site sequential distributed phosphorylation network ^ and recall the corresponding 
matrices 5("), y^"'\ £;("), H^ and M(") (cf. eqns. ^, (^, ^3^, [IJ^ and [153} ). The 
following are equivalent: 

(A) There exists a vector k G ^^q o-nd vectors a, b G -ff?>V with a ^ b satisfying hi 2a]) & \12b\) . 
that is 

S(") r(") {k, a) = 0, 5(") r^"' {k, b) = 0. 

(B) There exist vectors jj, G 12^"+^, /^ / 0, and {v, A) G IR^q x _K?JJ,, such that 



(C) There exist vectors fi G JR3"+3, ^ G -K^ and z^, A G 1R% x JR^'J,, mi/i ^ 7^ such that 
(a) the vectors fi and (^ satisfy 



y{n) ^^n(")e, 



(b) and the vectors v, A and ^ satisfy 



Xe]R%, free, 



and, for i = 1, . . ., n 



T^3i 



A3,e(2-)€i+(^-i)?^ 



Z^3i 



(16) 

(17a) 
(17b) 

(17c) 
(17d) 

(17e) 

(18) 

Concerning solutions a, b and n satisfying the polynomials (|12ap & (|12bp we observe the following 
(see also [S] Remark 7]): 

Remark 4.3 (Solutions a, b and k of the polynomial equations ()12ap fc (|12bp ). 
Recall the matrices M^'^' and E^"'' (cf. eq. (|15cp and eq. ()13bp ) and the monomial function $'"' 
defined in eq. (|lUal (cf. Remark I3.4p . Fix any two vectors /U G im (Af^"') and A G M^q. Then vectors 
(k, a) and [k, b) satisfying the polynomial equations 

5(")r(K,a) =0, 5(")r(K,6) = 0, 

are given by 

o, free in iR^^^, (19a) 

6 := diag (e^) a, (19b) 

K := diag ($(") (a-i)) ^(") A. (19c) 



i^3J-2 — A3j_2 T — , 

, Z^3i 

i^Si-l — A3j„i - — . 

('Dj There exists a vector /i G IR^^^^ , /i 7^ iwii/i 

;U G im iW^") 
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Note that, in particular, the vector a € IR^q is free and that a and b satisfy by construction 

In 6 - In a = /i G im (m^^'A . (20) 

This follows from the proof of Theorem 14. 2[ given in Section [5] below. □ 

As a consequence of Theorem 14.21 any two solutions a and b of the polynomial equations ()12ap 
& ()12bp satisfy the condition (|20p . From Definition 13.51 follows that for multistationarity a and 
b additionally have to satisfy Z'"' {b — a) = and hence b — a £ im(5^"'^). The existence of 
vectors a and b satisfying both conditions is covered by Corollary 14.51 given below (this result itself 
is a consequence of [4^, Lemma 1]; see also [3] and, for example, [lOj for an earlier, more informal 
discussion) . The corollary is based on the sign patterns of vectors and sets thereof. 

Definition 4.4 (Sign patterns defined by vectors and linear subspaces). 
Let V E ]W^. A vector a G {-1,0, 1}™ with 

(Tj = sign {vi) and hence a = sign (v) (21a) 

is called the sign pattern of v. A linear subspace V C iR™ then defines a set sign (V) of sign patterns 
in the following way: 

sign{V) ■.= {a e{-l,0,l}'^ \3v£V with sign (ct) = sign (u)} . (21b) 

Now we can state the result: 



Corollary 4.5 (cf. Lemma 1 of [3]). 

?3n+ ' 
''>0 



Two vectors a, b £ ^-^n^ with 



a^b (22a) 

ln6-lnaGim(M(")) (22b) 



6 - a G im (5(") j (22c) 

exist if and only if 

XW := sign (im (m^"))) fl sign (im (5("))) / {0} . (23) 

In case l[23\) holds, every pair /i G im (M'"'-'), s G im (S"'"-*) with 

sign (/x) = sign (s) / (24) 

defines a pair of positive vectors a and b with the property ()22p via: 

a = (ai)i=i,...,3n+3 (25a) 

with 



ef'i-l' 



if f'i^O 



(25b) 



l^aj > 0, free, 
6 = diag(e^)a. (25c) 

Proof. Apply l^, Lemma 1] with a as p, b as q, im (M'"-*) as Mi and im (5'"^) as M2. D 

As a consequence one can state a condition for multistationarity in system (|10bp defined by the 
network ([3|) based on the sign patterns defined by the linear subspaces im (5*^"^) and im (M^"-*) and 
the set Ad^""' . We will argue in Remark 14. 71 that this condition is satisfied if and only if at least one 
out of 2 (3^"^^ — 1) linear systems is feasible. Hence one may in fact establish multistationarity for 
the polynomial system (jlObp by an analysis of linear inequality systems. 
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Theorem 4.6 (A sign condition for multistationarity of systems (llObp ). 

Consider the n-site sequential distributed phosphorylation network ^ and recall the corresponding 
matrices S'("), Z(") and M("), (cf. eqns. (7^, (^) and H^) . There exists vectors a, b £ JR^"o+^ 
and K € IR^q satisfying 



5(") r(") {k, a) = 0, 5(") r^") (k, b) = 0, Z^") {b - a) = 0, 



7W(") = sign fim (m^"))) n sign fim (s'-'^A) / {0}. 



if and only if 



Proof Suppose that TW^") / {0}. Then, by Corohary [33] vectors a and b satisfying ()22ap - (|22cll 
are given by ()25ap - ()25cp . By Remark 14.31 (where a > is free and may be chosen as in (|25ap & 
(|25bp ) and Theorem 14.21 there exists a vector k € JR^*^ given by eq. (|19cp such that a, b and k are 
solutions of the polynomials (|12ap & ()12bp . By (j22cp a and 6 additionally satisfy Z*^"-* (6 — a) = 0. 
Vice versa, let a, 6 G JR^^q^^ with a / 6 and k G iR^To be given, such that (fT^ . (flj^ & (fT^ 
hold. From eq. (jlSp of Theorem 14.21 follows that ^ := ln6 — Ina G im (M*-"') (cf. also Remark 14. 3| 
eq. (I20D) and, as Z(") (6 - a) = ^ 6 - a G im (5(")), Corollary SJ implies that ([23]) holds. D 

Remark 4.7 (The condition TW^") / {0}). 

(A) 7V4'") 7^ {0} is equivalent to feasibility of (at least) one linear system: We want to emphasize the 
relation between condition ()23p and linear inequality systems: an element a G {—1,0,1}^"''''^ 
satisfies ()23p if and only if the following linear system is feasible: 

3s G JR3"+^ and S, e M^ such that 

Z(") s = 0, with ai Si > if 0-^ / and Sj = else, and (26) 

fi = M*^") ^, with (Tj ^j > if (Tj / and ^i = else. 

(B) Elements of { — 1, l}3"-+3. Recall the formulae (|25ap - ()25cp and observe that fn = <^ ai = bi. 
Hence a vector // (and corresponding s) with sign (^) G { — 1, lp"+3 yields a and b that differ in 
every coordinate, while a vector /i (and corresponding s) with sign (/i) G { — 1, 0, l}3"+3^ /U 7^ 
yields a and 6 that may be equal in some coordinates. 

(C) Testing 7W(") / {0}: Note that there are 33('"+^) (23("+i)) sign patterns a G {-l,0,l}3"+3 
(o" G { — 1, 1}^"^'^). A naive algorithm for computing elements a satisfying (|23|) would therefore 
be to enumerate ^ (3*-'^"'^^^ — l) (2''^'^+^'') of these sign patterns and check feasibility of every 
linear system (|26p associated to one of those. (As feasibility of ()26p for an element a implies 
feasibility for —a, one has to check only half of the sign patterns.) For the general problem 
of computing all sign patterns defined by two linear subspaces, a more advanced algorithm 
involving mixed integer linear programming is described in the thesis |29| . This software has 
aided the discovery of the sign patterns described below. □ 

4.2. Sign patterns for feasibility of system ()26p . First we turn to the subproblem Z^'^> s = 
with sign (s) = cj^"^' in ()26p for a given cj^""' G { — 1,1}^"'^^. To this end, we consider sign 
patterns in column and matrix form. With (n + 1) column vectors T(n\ = (tio,T2o,T3o) , •••, 
''"(n) = {Tin,T2n,Tsn)'^ from {—1, 1}^ we associate the sign pattern matrix 

r=(r(o)---r(„))G{-l,l}3x("+i). (27) 

and the a sign pattern column t^^' = coI{t) = col{TtQ\, ...,T(„\) G { — 1, l}3"+3. 
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Lemma 4.8 (Feasibility of Z(") s = and sign (s) = o-(") G {-1, 1}^"+^). 
For a given sign pattern matrix 

(o"io ••• crin \ 
O"30 ... 0-3n / 

and its associated sign pattern vector o"'"-* = col{crrQ\, ...,ar„\) S { — 1,1}'^""^^, the subproblem 

Z(") s = with sign (s) = cj(") (28) 

m \20\) is not feasible if and only if a has one of the following three properties: 

(1) The elements of the first row are all of the same sign cJio- 

(2) The elements of the second row are all of the same sign (T20 with (72oO"io < and ct2oO"3o < 0. 

(3) The elements of the third row are all of the same sign a^Q. 

Proof. The claim follows easily from the fact that Z^") s = is equivalent to the decoupled linear 
system 

n n n 

y^ Slj =0, ^ S3j =0, ^ S2j = Sio + S30 . 
j=Q j=0 j=0 

D 

We now turn to those sign patterns o"'"^ E { — 1, l}3"+3^ where (|28p is feasible and where, additionally, 
the second subproblem 

/i = m(") e with sign (/u) = fj^") (29) 

is feasible. In particular, we will show that (|28p and (|29p are only feasible for sign patterns involving 

C7i:=(+1,+1,+1)^, a2:={+l,-l,+lf, ^3 := (-1, -1, -1)^ = -ai (30a) 

a4 :=(+!, -1,-1)^, a5:={-l,+l,-lf = -a2, ag := (-1, +1, +1)^ = -^4, (30b) 
while it is infeasible for sign patterns involving 

aj := (+1, +1, -1)^ , as := (-1, -1, +1)^ = -aj. (30c) 

We use these aj to encode sign patterns a^"^' G { — 1, l}3"+3 Pqj- example, for n = 2 we encode 

a("):=( 1, 1, ![ 1,-1, l\-l,-l,-lf eM^ 

as the 3x3 array 

a = (cri(T2 0-3) 
with the columns ui, cr2 and (T3. We use the symbol [r]j as shorthand notation for the i-fold 
concatenation of r E { — 1, 1}^, i.e., for the 3 x i-matrix 

[A = {r_^). (31) 

i times 

For an example with n = 5, 

^ = (kl]3 [0-2]2 0-3) = (0-i(Ti0-iCT2CT20-3) 

stands for the 18- dimensional sign pattern 

a(") = ( 1, 1, 11 1, 1, 11 1, 1, 11 1,-1, 11 1,-1, 11-1,-1,-1)^ 



"-5 -^5 



In the following Theorem 14.91 we characterize the sign patterns a^"'' € {—1, l}3"+3 tj^^t are feasible 
for (j26p for n > 2. In Section [6] below we will prove this result. 
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Theorem 4.9. Recall the linear system i26\} associated to the sign pattern vector a^"'' € {—1, 1}^"+^ 
and hence to the equivalent sign pattern matrix a G { — 1, 1}3>^("+1), Por n > 2 the linear systems 
i26]} associated to a^'^' and to — o"'"' are feasible (and hence Ai^'"'' from (I23p is non-empty) if and 
only if a is one of the following sign patterns (jsij) to (jsyj); 

• For n > 2: 

a = ((T2 [cTsJi-^ kilja) ^*^^ n, ^2 > 1 and ii + 12 = n (si) 

(^ = (0-3 kslii kilia) ^^^^ n, «2 > 1 and ii-iri2 = n (32) 

a = (0-4 [0-5]^^ kilia) ^^^^ n = I7 -^2 > 1 and i^ + ig = n (S3) 

a = (0-4 [0-3] .^ kilja) ^*^^ ^i> ^2 > 1 and ii + 12 = n . (34) 

• Additionally for n > 3: 

a = (0-2 [0-3] -^ [cr2]j2 kijig) w'ii/i ii, 12, i3 > 1 and i^ + ig + «3 = "- (ss) 

o- = (0-3 kalii kslia kilig) w'*^^ n, «2, is > 1 an^c? i^ + ig + «3 = "- (se) 

0- = (0-4 [0-3] -^ [0-5]^^ kilig) with ii , i3 > 1; «2 = 1 and i^ + ia + 13 = n . {37) 

As an immediate consequence of Theorem 14.91 above we note the following: 

Remark 4.10 (Multistationarity in n-site sequential distributive phosphorylation for n > 2). 
As a consequence of Theorem 14.61 and 14.91 we conclude that for every n > 2 there exists a vector 
of rate constants k such that the system ()12ap . (|12bp & ()12cp has at least two distinct positive 
solutions a and b. Every element a described in Theorem 14.91 defines a feasible linear system 
of the form (|26p . And solutions thereof parameterize the aforementioned rate constants via the 
equations ()25ap . (|25bp and (I19cp . It will be interesting to systematically explore the corresponding 
regions in the rate constant space, especially with respect to biological plausibility and function. □ 

Example 4.11 (The sign pattern (JS2J))- 

For 2 < n < 5 one obtains from the formula ([S2J) the following patterns: 

n = 2 n = 3 n = A n = 5 

0-3 <75 CTs 0-5 (T5 (Ti 

_,.-—- 0-3 (T5 (T5 (75 (Ti ___^ 

__— — 0-3 0-5 (75 CTi ____^ 0-3 0-5 (T5 (T5 (Ti Ui 

^■3 0-5 <7l -_ ____—- (73 (T5 (T5 (71 CTi- ___^ 

^~-- (T3 (T5 (Tl (Ji C:;! ___lr= 0-3 (T5 (75 (Tl (Tl (Tl 

^-— (T3 0-5 (Tl (Tl fTi. ____^ 

(T3 (T5 (T5 (Tl (Tl (Tl 

We summarize for ^s^ in case of n = 5: 

• The integer 5 can be partitioned in the sum of two integers ii, 12 in the following two ways: 
5 = 1+4 = 3 + 2. As we have to take the ordering of ii and ^2 into account (and hence, for 
example, (4, 1) differs from (1,4)) we obtain the following pairs of integers: 

{(1,4), (4,1), (3,2), (2,3)}. 

• This yields the sign patterns 

0-3 k5]4 kl]l> 0"3 ksli [o-l]4. Cr3 [a^]^ [cJi]2 , ^3 [a^]2 [CTllg 

□ 
Concerning the number of elements a^^' € { — 1, l}3"+3 corresponding to (jsij) - (jsyj) we observe 

Proposition 4.12. 

For n > 3 there exist (n — l)(n + 2) elements in { — 1, l}3x("-+i) ^^ ^/^g ^orm (^1]) - Jst]) . 
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Proof. We note that for n > 3 fixed formulae ( jsg] ) and ( jggj ) each yield as many elements as there 
are partitions of n into three positive integers (taking order into account). Similarly (pi]), ^s^ and 
( [S4] ) each yield as many elements as there are partitions of n into two positive integers (taking order 
into account), while ( [57] ) yields as many elements as there are partitions of n — 1 into two positive 
integers (taking order into account). Finally, ([S3J) yields one element. 

First we count the number of ways an integer can be partitioned into the sum of two positive integers 
(if order is taken into account). One has 



n 



(n-l) + l 



[n 



2) + 2 



(n- (n- 1)) + (n- 1). 



Note that there is no other way to partition n into the sum of two positive integers. Hence we 

conclude that there is a total of n — 1 ways to partition n into two positive integers. 

Second we count the number of partitions of n into three positive integers. To this end we observe: 



n 



(n-2) + (l + l) 

(n — 3) + {2 ways to write 3 as the sum of two positive integers} 



= {n — j) + {j — 1 ways to write j as the sum of two positive integers} 

= (n — (n — 1)) + {n — 2 ways to write n — 1 as the sum of two positive integers} 
Consequently there is a total of 

n-l 



Y,U-l) = -{n-2){n-l) 



i=2 

ways to partition n into the sum of three positive integers. Now we can sum over the contribution 
of each of the formulae 

(n-2)(n-l)+ 3 (n-l) + (n - 2) + ( 1 ) = (n - l)(n + 2) . 



( tsst & imj in}' BiJ ^ 61J inJ 



t^ 



D 



5. Proof of Theorem 14.21 

To prove the Theorem 14.21 we show the following implications: 

(A) ^ {D) ^ iO => (B) ^ (A) 

(A) ^ jD) : 

Suppose there exists a, b € M 

Theorem 4.3 of 1201 one has 



i3n4 
>0 



and K G M^^Q^^ 



In 6 



such that (fT2al) & (IT2b]l hold. Then, by 
In a G im [A'^) , 



with 



A 



1 



1 




1 1 1 


2 2 2 




1 1 


1 1 




1 1 1 


1 1 1 





n n 




1 



n 
1 
1 



where the columns of A have been reordered to match the ordering introduced in Table [T] (in the 
aforementioned reference ^20j . the columns of A correspond to the exponents of ti, t2, t^ in the 
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statement of Theorem 4.3, the matrix itself is displayed in the proof of said Theorem). We want to 
show that im (M^"') = im [A^^ . To this end we read off the following basis of ker (A): 

(n) (n) (n) , (n.) (n) (ra) , (n) (n) . -, 

which is also a basis of ker ( M^"' J. As both subspaces im (M*-"-*) and im (^"^) have the same 
dimension we therefore conclude 

im ('m^")) = im {A^)^ and hence im [m^")) = im (A), 

implying (A) =^ (D), with ^ := In 6 — In a. 

(^) ^ (g) : 

We first show that 3^(")^ M(") = n("). To this end recall the definition of 3^(")^ from eq.([8d| 



(n)J 



>b(lf 06.1x3 


06-2x3 


06-3x3 






06-nx3 


^o(2)" 


>^o(3)^ 










Y^(n-\Y 








y^(nf 









y 



Now it is straightforward to compute that 

r M„ (0) " 

Y^iif : =[0 (-i + 2)l {i-\)\\=:A{i), 

_ Mnii) _ 

where, as in Definition 14.11 we use the symbol \ to denote a 6- vector filled with the number 1 and, 
accordingly, we use the symbol to denote a 6-vector filled with the number 0. Thus one arrives at 

A{\) 



yi.nY jVfW 



A(n) 



Further note that 



A(S) Uij =[H + 2)1 (^-l)l] (|)=no 



(i) 



6 



(cf. Definition 14.11 and in particular equation ()14ap ). Thus, given any nonzero ^ G im (M^"-*) as 



v6y 



one arrives at 



yW M^""^ \ix 



no(i) 



noN 



n(") i. 



Concerning v and A as in (|17bp &: (I17cp - ()17ep we note that ^i, ^2 from above together with any 
positive A can be used to obtain v as in ()17cp - ()17ep . Hence we conclude {D) =^ (C). 

(g) ^ (i?) : 

Let ^ G JR3"+3 with /i 7^ and ^ G JR^ be given, such that (|17ap holds. Further let v and A be 
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given that satisfy (flTbl) & (fTTcll - (fTTe]) . that is A G JR^^g and 

i^3i = A3ie(^"*)^i+(*"^)^^ z^3j_i = Asj-i -^ and z^3j_2 = A3j_2 -^,i = 1, . . . ,n. 

Observe that these i', A satisfy 

, I^3J-2 , l^3i-l 1 Z^3i /„ -N /- I r 1^ t 
In = In = hi -— = (2 - t) 6 + (« - 1) 6 

A3i-2 A3j_i A3j 



and hence 

Define 
and let 



m — = in - — and in — = in - — . 

A3i-2 + A3j A3J Asi-i + Asi Asj 

^(i) '■= (^3i-2, J^3i-1, ^3i) and A(j) := (A3j_2, A3i_i, Asi) 



i^ = col (i/(i), . . . , !/(„)) and col (A(i), . . . , A(„)) . 
Recall the matrices E^^' and E from eqns. (|13bp & (|13ap and observe that 



E(»)A V E\i, EA,„, 

and 

1 -^ ^(*) A ^3i~2 + J^3i 1 Z^3i-2 1 i^3i , Z^3i-1 + Z^3i , '^3i-l , Z^3i 

-^ ^(j) V A3j_2 + Asi A3J2 A3j A3j_i + A3j A3j_i A3J 

Hence one has for z/. A from above 



,|M.[p_.u(.-i,i](|) 



Hi) 

and 

S(") A 

Consequently, if there exist vectors fi € iR'^"'*'^, /i 7^ and ^ € iR^ satisfying ([TTaj) and vectors z/, A 
satisfying ()17bp . ()17cp - ()17ep . then /x, z^ and A also satisfy ()16p . Note that z^ and A satisfying ()17bp 
& ()17cp - (I17ep are positive. Thus we conclude (C) =^ (i?). 

JB) ^ (^) : 

Finally, assume /x G iR^n+s ^ith ^ / and (z^, A) G iR^To x IR^^ satisfy ([16]). Fix a G iRI^g and 

define 

b := diag (e'") a 

K := diag ($(") (a-i)) ^("^ A. 

Then 

5(") diag (k) $(") (a) = 5(") diag ($(") (a-^)) diag (^("^ a) $(") (a) = 5^") .E^") A = 

and(as,by(Il6D, y")^^ = ln|^) 

5(") diag (k) $(") (6) = 5(") diag (k) diag ($(") (a)) e^*"'^/^ = 5^") .E^") z^ = 
and therefore (i?) =^ (A) where we have used equation (jllap of Remark 13.41 □ 
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6. Proof of Theorem 14.91 

First, we observe that sign (M^") ^) = cr^") implies sign (M^") (-^)) = -a^") and that Z^") s = 
implies Z^"^' (— s) = 0. Hence, if the system (|26p associated to o"'"^ is feasible for (^, s), then the 
system ()26p associated to — cr^"^ is feasible for (— ^, — s). Thus cr G A4^^' implies —a € Al^"'^ Hence 
the linear system (|26p is feasible for elements a^^' € { — 1, lp'^+3 g^^^ hence — cj'"-' if and only if a^"'' 
corresponds to one of the sign pattern matrices in ( jsij ) - (jsTJ) . 

We now turn to the necessity of the sign patterns (jsij) - ( jgTJ ) . To this end we parametrize the range 
of the regular matrix M (0, n) from eq. (J15ap by 



w{0) :-- 



Consequently, one has for M (i, 



y 




X 


= 


_y - z_ 





n) 



1. 



.71, 



-n + 1 n 

n —n 

-n + 2 n — 1 

from eq. (|15bp 



X + y 

X + y + z 



(32) 



■w U 



wi(z) 
tt;2(i) 



X + iz + y — z 

X + iz 
X + iz + y — z 



-i + 2 i-1 

n — i —n + i 

-i + 2 i-1 



X + nz 

x + y 

X + y + z 



(33) 



Observe that the components W2{i) and wi{i) = w^{i) = W2{i) + {y — z) are affine functions of i 
so sign changes can be easily read off. We now derive the necessary conditions for sign(t(;(0)) to 
be one of the four sign patterns cri , o"2 , (T4 and cry from (|30ap - ()30cp under the side condition that 
system ()28p is not unfeasible (cf. Lemma [4. 8p . The indices i^ and jk to follow will always be > 1. 

(1) sign (u;(0)) = £72, i.e., y > 0,x < 0,y - z > 0: 

For z < 0, the W2i are all negative so that one has o"2o = 0'2i = • • • = (T2n = —1 and 
fio = fso = 1- By item (2) of Lemma [4.81 z <0 cannot lead to a feasible sign pattern. For 
z > the possible sign pattern matrices are 

(0"2, kslii, [cnlia) ' {^2, kslji, [0-2]j2> Kljs) 

for ii+i2 = n and ii+j2+i3 = n where the cJi entries are necessary by item (2) of Lemma l4.8l 
Observe that the sign pattern (o"2, [cijn) is unfeasible by item (1) of Lemma [ 

(2) sign {w{0)) = 0-3, i.e., y < 0,x < 0,y - z < 0: 
For z < 0, the components of all the w{i) are of the same sign. By item (1) of Lemma [ 
z < cannot lead to a feasible sign pattern. For z > 0, the w^i) may generate the sign 
pattern matrices 



(0-3, kslii, [cnlia). (o"3. k3]ii, [o-5]j2, Klja) 

with ii + ^2 = n and ji + J2 + js = n where the cri entries are necessary (otherwise the first 
rows do not offer a sign change). The direct passage from (T3 to o"i is obviously impossible 
since W2{i) < implies wi{i + I) = W2{i) + y < W2{i) < 0. 

(3) sign(t(;(0)) = 0-4, i.e., x < 0,0 < y < z: 

At first, possible sign pattern matrices are 

for ii+i2 = n and J1+J2+J3 = ''^ where the o"i entries are necessary by item (3) of Lemma l4.81 
Observe that the sign pattern ((74, [cri]„) is unfeasible by item (1) of Lemma [4.81 Note that 
(T5 can appear just once because of wi^i^i = W2i + y > W2i- 

(4) sign(tt;(0)) = aj, i.e., x > 0,0 < y < z: 

Because of 2; > 0, the components wu are positive for i = 0, ...,n. By item (1) of Lemma[ 
(Ty cannot lead to a feasible sign pattern. 
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We finally establish the sufficiency of the sign patterns ( [si] ) - ( [S7J ) by proving that they are realizable. 
To this end, we present for each (sfc), k € {1, ..., 7}, a vector 

^ = col{x,x + y,x + y + z) e M^ (34a) 

leading to the sign pattern (sk) for M^")^: 



For (si) 
For (S2) 
For (S3) 
For (S4) 
For (ss) 
For (se) 
For (S7) 



x = i2-j, y = -n + \, z = l. 

X = n — J, y = — n — ii + |, z = l. 

X = n — J, y = —n+j, z=l. 

x = i2-\, y = -n + \, z = \. (34b) 

2; = Z3-|, y = -n + i2 + |, z = \. 

x = i2^h-\-, y = -n-i2 + j, z = l. 

x = i2 + l, y = -n + \, z = l. 

7. Discussion 



The main topic of this contribution is the existence of multistationarity for n-site sequential distribu- 
tive phosphorylation. Via the equations (|19c|) and (|25ap - (|25c|) we were able to link the existence of 
pairs of steady states (a, b) and a corresponding vector of rate constants k to solutions of the linear 
systems (|26p . These systems are uniquely defined by sign patterns a^"'' G { — 1,0, l}3"+3 a,nd hence 
the existence of a single sign pattern defining a feasible system (j26p is sufficient. Concerning sign 
patterns a G { — 1, l}3"+3 where (j26p is feasible. Theorem 14.91 establishes four such sign patterns for 
n = 2 and (n — 1) (n + 2) for n > 3. 

Unfortunately the biological interpretation of elements a where (j26p is feasible is unclear (unlike in 
the case of multistationarity itself, that is usually associated with a specific biological function). If 
different elements could be associated to different biological functions, this would limit the number of 
(different) biological functions a network can perform and hence be of significant biological interest. 
Moreover, knowledge of all sign patterns a defining feasible (j26p would then correspond to knowledge 
of all biological functions a network is capable of. 

However, the following observations concerning pairs of steady states (a, b) can easily be made 
by inspection of the feasible sign patterns ( pij ) - ([S7J) and the structure of the matrices M^*' from 
()15bp : (i) every sign pattern ends with the triplet (1,1,1) or ( — 1,— 1,— 1) and hence the last 
three components of steady state b are greater (smaller) than those of a (recall that ^ = In ^ and 
hence /ij > implies 6j > Oj and ^j < implies 6j < Oj). (ii) The first and the third row of the 
matrices M^"^' are identical. Hence fi^i+i = /^si+s, i = 1, . . . , n. For a pair of steady states (a, 6) 
these correspond to the ratio of steady state concentrations of kinase-substrate complex Ai-ipEi (in 
case of ^3j+i) and phosphatase-substrate complex Ai^ E2 (in case of /U3J+3, cf. Table [1]). Hence, for 
any pair of steady states, the ratio of the steady state concentrations of kinase-substrate complexes 
equals that of phosphatase-substrate complexes. 

Finally we'd like to discuss two observations related to the formulae (|19cp and (|25ap - (|25cp : consider 
a sign pattern a as in ( jsij ) - ( [57] ) and pick vectors /i G M^^\ s G S^^' with sign (s) = sign (^). Let 
{a^ , b^) be the pair of steady states defined by //, s via (|25ap - (|25cp . Let {a~^ , b~") be the pair of 
steady states defined by —jj, and —s. Then it is not difficult to see that a~" = b'^ and b~'^ = a" (cf. 
a similar discussion in [1]). 

Further note that a given pair (a, b) obtained via (|25ap - (j25cp is a pair of steady states for all 
vectors of rate constants defined by the formula (|19cp . where A G JR^q, free. This may be interesting 
for two reasons: first, it allows to replace the 6n rate constants /cj, li by the 2>n 'coordinates' Aj 
and hence a reduction of the degrees of freedom in the dynamical system (jlObp . Second, it means 
that multistationarity is robust with respect to perturbations within the image of the map (|19cp 
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and fragile with respect to transversal perturbations. In a future publication we will explore the 
image of this map with the aim of establishing parameter values for multistationarity in biologically 
meaningful domains. 
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